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Abstract 

The Newman-Penrose formalism is used to deal with the quasinormal modes(QNM's) 
of Rarita-Schwinger perturbations outside a Reissner-Nordstrom black hole. We ob- 
tain four kinds of possible expressions of effective potentials, which are proved to 
be of the same spectra of quasinormal mode frequencies. The quasinormal mode 
frequencies evaluated by the WKB potential approximation show that, similar to 
those for Dirac perturbations, the real parts of the frequencies increase with the 
charge Q and decrease with the mode number n, while the dampings almost keep 
unchanged as the charge increases. 
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For several decades, the QNM's of black holes have been of great interest 
both to gravitation theorists and to gravitational- wave experimentalists [1,2,3] 
due to the remarkable fact that QNM's allow us not only to test the stabil- 
ity of the event horizon against small perturbations, but also to probe the 
parameters of black hole, such as its mass, electric charge, and angular mo- 
mentum. QNM's are induced by the external perturbations. For instance, if 
an unfortunate astronaut fall into a black hole, the surrounding geometry will 
undergo damped oscillations. They can be accurately described in terms of a 
set of discrete spectrum of complex frequencies, whose real parts determine 
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the oscillation frequency, and whose imaginary parts determine the damped 
rate. Mathematically, they are defined as solutions of the perturbation equa- 
tions belonging to certain complex characteristic frequencies which satisfy the 
boundary conditions appropriate for purely ingoing waves at the event horizon 
and purely outgoing waves at infinity [4]. 

Recent observational results suggest that our Universe in large scale is 
described by Einstein equations with a cosmology constant. Motivated by 
the recent anti-de Sitter (AdS) conformal field theory (CFT) correspondence 
conjecture[5], much attention has been paid to the QNM's in AdS spacetimes[2,3] 
The latest studies on quantum gravity show that QNM's also play an impor- 
tant role in this realm due to their close relations to the Barbero-Immirzi 
parameter, a factor introduced by hand in order that loop quantum gravity 
reproduces correctly entropy of the black hole[6,7,8,9]. 

A well-known fact is that quasinormal mode (QN) frequencies are closely 
related to the spin of the exterior perturbation fields [10]. Previous works on 
QNM's in black holes has concentrated on the perturbation for scalar, neu- 
trino, electromagnetic and gravitational fields[3]. However, few has been done 
for the case of Rarita-Schwinger fields, which closely relate to supergravity. 
According to the supergravity theory, the Rarita-Schwinger field acts as a 
source of torsion and curvature, the supergravity field equations reduce to 
Einstein vacuum field equations when Rarita-Schwinger field vanishes. This 
can be seen from the action of supergravity, namely [11], 



where L2 and L3/2 represent the Lagrangian of gravitational and Rarita- 
Schwinger fields, respectively. We hence expect to obtain some interesting 
and new physics by investigating ONM's of Rarita-Schwinger field. 

As to QNM's, the first step is to obtain the one- dimensional radial wave- 
equation. We usually start with linearized perturbation equations. Two often 
used ways are available for obtaining the linearized perturbation equations. 
One is a straightforward but usually complicated way that linearize the Rarita- 
Schwinger equation directly and deduce a set of partial differential equations 
(one can see Ref.[12] for details); The other is provided by the Newman- 
Penrose (N-P) formalism, which end up with partial differential equations in 
r and 6 instead of ordinary differential equations in r. Torres [13] have deduced 
the linearized Rarita-Schwinger equation in N-P formalism for a type D vac- 
uum background. We start with the Rarita-Schwinger equation in a curved 
background space-time 



or, in the Newman- Penrose notation, namely, the Teukolsky's master equations [14, 15] 
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{[D - 2e + e* - 3p - p*] (A - 87 + //) 
-[5 - 2/3 - a* - 3r + 7r*](^ - 3q; + tt) - *2}^3/2 = 0, 



(2) 



and 



{[A + 27 - 7* + + p*]{D + 3e - p) 
-[5 + 2a + /5* + 37r - r*](5 + 3/3 - r) - ^2}$-3/2 = 0. (3) 

Here we have introduced a null tetrad {l'^ , n'^ , m'^ , fh'^) which satisfies the or- 
thogonality relations Ij^n^ = —m^fh^ = 1 and /^m^ = /^m^ = n^m^ ~ 
n^fh^ — 0, and the metric conditions g^^, = l^ny + n^l^, — m^fhi, — fh^m^. Ac- 
cording to these conditions, we can take the null tetrad as = (e-^^^^), 1, 0, 0), 
^ 1(1, -e^^W, 0, 0), = ^ (o, 0, 1, m'^ = ^ (o, 0, 1, . The 

corresponding covariant null tetrad is = (1, — e^^'^^''^ 0, 0), n^^ = \{e'^^^'^\ 1, 0, 0) 
, = — (0, 0, 1, zsin6'), fh^ — — (0, 0, 1, — isin^). The non-vanishing 
spin coefficients read 

^ r 2V2r 2r ' 4^ ^ ' ^ ^ 

and only one of Weyl tensors is not zero, i.e., ^2 = — ^^'^ ' where a prime 

denotes the partial differential with respect to r. 

In standard coordinates, the line element for the Reissner-Nordstrom space- 
time can be expressed as 

ds-" = -e^^^'Ue + e-^^^'-^dr^ + (dO'' + sin'^ed^'') , (5) 

with 

e-"' = l-^ + ^, (6) 

where M and Q are the mass and charge of the black hole, respectively. 

The directional derivatives given in Eqs.(2) and (3), when applied as deriva- 
tives to the functions with a t- and a (/^-dependence specified in the form 
^t{u}t+m(p) ^ become the derivative operators 



D = %, A = -A^t^ '^ = ^^o^ '^=^^^0, (7) 



where 



^ „ iur'^ A' X ^ iur'^ A' 

^n = dr + +n - —, S>1^ dr- ^ + n- —, 

^n^de + —-: + ncote, ^ de - —- + ncotO, (8) 

sm 6 sm 
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and 

A = r^-2Mr + Q^ (9) 

It is obvious that ^„ and ^| are purely radial operators, while =Sf„ and 
are purely angular operators. After some transformations are made, Eqs.(2) 
and (3) can be decoupled as the two pairs of equations[10]. 



[A^_i/2^o^-4^t^r]P+3/2 = 

-^-1/2-^3/2^+3/2 = 



-^-^+3/2, 



(10) 

(11) 



and 



[A# ,/2^0 + W]P_3/2 = AP_3/2, (12) 
-Sf-1/2-^3/2^-3/2 = -AA_3/2, (13) 

where A is a separation constant. The reason we have not distinguished the 
separation constants in Eqs.(ll)- (13) is that A is a parameter that is to be 

determined by the fact that A_|_3/2 should be regular at ^ = and 6 = n, and 
thus the operator acting on v4_3/2 on the left-hand side of Eq.(13) is the same 
as the one on A^2/2 in Eq.(ll) if we replace 6* by tt — ^. 

In Reissner- Nordstrom black hole, the separation constant can be determined 
analytically [1 5 ,16,17] 

f(/ + 3)(/ + l) for j = l + s, 

\l{l-2) for j = ^ ' 

where / = 2,3,4,---. Note we only consider the case for j — I + s in our 
following discussions, the case for j = Z — s can be easily obtained in the same 
way. Since P+3/2 and P-3/2 satisfy complex-conjugate equations (10) and (12), 
it will suffice to consider the equation (10) only. 

2 

By introducing a tortoise coordinate transformation dr^ = ^dr, and defining 
A± = ^ ± ia;, F = r~^P+3/2, one can rewrite Eq.(lO) in a simplified form 

AV + PA_F - gy = 0, (15) 

where 

d r-6 _ A r 9A A' 1 

(16) 



Transformation theory [4] shows that one can transform Eq.(15) to a one- 
dimensional wave-equation of the form A'^Z = VZ by introducing some param- 
eters (certain functions of r^, to be determined) (^(r*), x(^*)) /?i(^*)) ^i(^*)i and 
several constants (to be specified) P2, T2, k, Ki- If we assume that Y is related 
to Z in the manner Y = ^VZ+TA+Z, and the relations T ^ Ti{n) +2iuj, P ^ 
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f3i{r*) + 2iujP2, one can obtain a equation governing /32, k, ki derived from 
eq.(15)[10] 



^3/2 



F-P2 



+ 



(17) 



where we have defined F = -^21 and ', r,,,' denotes the differential with respect 
to r*. A key step to obtain the expression of potential V is to seek available 
/32, and ki whose values satisfy Eq.(17). Further study shows that F as 
defined in the previous text does satisfy Eq.(17) with the choice 



/?2 = ±2Q, K = 0, Ki = ±A\/1 + A. 
The potential V can then be expressed as 



V ^ 



A3/2 



Kl){KiF - jS^F^r.) 



{F-mF' 



HI) 



(18) 



(19) 



Note that the sign of (i^ and k\ can be assigned independently, so there exists 
four sets of expressions of V . Here we denote them by V\3 = 1, 2, 3, 4). Chan- 
drasekhar pointed out that there has a closed relation between the solutions 
belonging to the different potentials showed in equation(19) (see Ref.[4] §97(d) 
for details) 



where 



F + P2 



F + P? 



lUJ- 



Pi 



k{F - piF' 



F + Pi 



F2 



PI 



dr* 



(20) 



p\ 



[F-mF-m 

One can easily prove that the potentials vanish when we let r 



±00 



g2M 
„-2 



-00, 



+00. 



(21) 
(22) 



A direct consequence of this property is that the wave-function has an asymp- 
totically ffat behavior for r — > ±00, i.e., Z — e^*'^'* (this is just the boundary 
conditions of QNM's). It has shown that in asymptotically flat spacetimes, so- 
lutions related in the way showed in Eq.(20) yield the same reflexion and trans- 
mission coefficients, and hence possess the same spectra of QN frequencies [4]. 
Moreover, we can easily obtain the potentials for negative charge by rearrang- 
ing the order of for positive ones since j32 equals to ±2(5. Therefore, we shall 
concentrate just on potential with /?2 = '2Q{Q > 0, say) and Ki — Xy/1 + A in 
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our following works. 

We hence know that the radial equation (10) can be simplified to a one- 
dimensional wave-equation of the form 



dr'i 



+ u?Z = VZ, 



(23) 



where 



V ^ 



(F-2Q)(F2-4Q2 



(24) 



Note that we have written as V because we only work with one case of the 
potentials. 

The effective potential V{r, Q, I), which depends only on the value of r for 
fixed Q and I, has a maximum over r e (r+,-|-oo). The location tq of the 
maximum has to be evaluated numerically. An interesting phenomenon is that 
the position of the potential peak approaches a critical value when / — > oo, 
i.e., 

ro{l ^oo)^ 3M. (25) 

Obviously, the effective potential relates to the electric charge of black hole. 
Figure 1 demonstrates the variation of the effective potential V{r, Q, I) with 
respect to charge Q for fixed 1 — 2. Prom this we can see that the peak value 
of the effective potential V increases with Q, but the location of the peak 
decreases with charge. This is quite consistent with the case for Dirac pertur- 
bation in Reissner-Nordstrom black hole spacetimes. 

We now evaluate their frequencies by using third-order WKB potential 
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Fig. 1: Variation of the effec- 
tive potential for massless grav- 
itino with M = 1, Z = 2, Q = 
0,0.3,0.6,0.9 



approximation [10], a numerical method devised by Schutz and Will[18], and 

was extended to higher orders in[19,20]. Due to its considerable accuracy for 
lower-lying modes[21], this analytic method has been used widely in evaluating 
QN frequencies of black holes. Noting that during our evaluating procedures, 
we have let the mass M of the black hole as a unit of mass so as to simplify 
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Tabic 1 

QN frequencies of Rarita-Schwinger field in RN black hole for I = 5 (third-order 
WKB approximation) 



Q 




n = 




n = ' 


1 




n = 2 




n = 3 




n = 4 





1, 


.3273+0.09581 


1 


.3196+0, 


,28811 


1 


.3048+0.48241 


1 


.2839+0.67951 


1 


.2582+0.87941 


0.1 


1 


.3295+0.09591 


1 


.3218+0, 


.28831 


1 


.3070+0.48271 


1 


.2863+0.67981 


1 


.2606+0.87991 


0.2 


1 


.3364+0.09601 


1 


.3287+0, 


.28881 


1 


.3140+0.48351 


1 


.2933+0.68091 


1 


.2678+0.88131 


0.3 


1 


.3481+0.09631 


1 


.3405+0, 


.28951 


1 


.3260+0.48471 


1 


.3055+0.68261 


1 


.2802+0.88341 


0.4 


1 


.3653+0.09661 


1 


.3579+0, 


.29061 


1 


.3436+0.48641 


1 


.3234+0.68491 


1 


.2985+0.88631 


0.5 


1 


.3890+0.09701 


1 


.3817+0, 


.29181 


1 


.3677+0.48841 


1 


.3481+0.68761 


1 


.3238+0.88961 


0.6 


1 


.4206+0.09741 


1 


.4136+0, 


.29301 


1 


.4001+0.49031 


1 


.3811+0.69021 


1 


.3576+0.89281 


0.7 


1 


.4627+0.09771 


1 


.4560+0, 


.29381 


1 


.4432+0.49161 


1 


.4251+0.69171 


1 


.4027+0.89451 


0.8 


1 


.5197+0.09761 


1 


.5135+0, 


.29321 


1 


.5016+0.49041 


1 


.4848+0.68981 


1 


.4639+0.89161 



the calculation. The values are hsted in Table 1, where we only list the values 
for / = 5 as an example. Values for other mode numbers can easily obtained 
in the same way. As a reference, we have also evaluated the values (listed in 
Table 2) for / = 5 by using the first-order WKB potential approximation [18]. 
Obviously, compared to the first-order approximation, great improvement, es- 
pecially for larger n, has been made for third-order approximation. 

Figure 2 demonstrates the variation of real and imaginary part of the QN 
frequencies with different Q and n for / = 5. It shows that the real part of 
the quasinormal mode frequencies increases with the charge Q, while decreases 
with n. But things are totally different for the imaginary part as showed in the 
figure, whose values almost keep unchanged as the charge increasing, whereas 
them increase very quickly with the mode number. Furthermore, there is also 
an interesting phenomena that the larger the charge is, the smaller effect of n 
on the real part of QN frequencies may have. 
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Tabic 2 



QN frequencies of Rarita-Schwinger field in RN black hole for I = 5 (first-order 
WKB approximation) 
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28321 


1 


4282+0 45691 


1 


4862+0 61471 


1 


5507+0 


75751 


0.4 


1 


3737+0.09661 


1 


3995+0 


28431 


1 


4451+0.45891 


1 


5031+0.61771 


1 


5677+0 


76151 


0.5 


1 


3973+0.09701 


1 


4229+0 


28571 


1 


4683+0.46141 


1 


5262+0.62151 


1 


5908+0 


76661 


0.6 


1 


4288+0.09741 


1 


4541+0 


28711 


1 


4991+0.46411 


1 


5567+0.62571 


1 


6213+0 


77241 


0.7 


1 


4707+0.09771 


1 


4954+0 


28821 


1 


5398+0.46641 


1 


5968+0.62971 


1 


6610+0 


77831 
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1 


5273+0.09751 


1 


5512+0 


28801 


1 


5942+0.46711 


1 


6500+0.63181 
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7133+0 


78231 
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Fig. 2. Variation of the QN frequencies for V with different Q and n for Z = 5. 
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